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Abstract. 

In this paper we introduce a generalized concept of field-transforming 
metamaterials, which perform field transformations defined as linear relations between 
the original and transformed fields. These artificial media change the fields in 
a prescribed fashion in the volume occupied by the medium. We show what 
electromagnetic properties of transforming medium are required. The coefficients of 
these linear functions can be arbitrary scalar functions of position and frequency, which 
makes the approach quite general and opens a possibility to realize various unusual 
devices. 
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1. Introduction 

It has been recently realized that metamaterials - artificial electromagnetic materials 
with engineered properties - can be designed to control electromagnetic fields in rather 
general ways. The concept of "transformation optics", which is based on finding 
artificial materials that create the desired configuration of electromagnetic fields, has 
been developed by several research teams, see e.g., [U El [3] 

In the known approaches one starts from a certain transformation of spatial 
coordinates and possibly also time, which corresponds to a desired transformation of 
electromagnetic fields. It has been shown that transformation of spatial coordinates 
(accompanied by the corresponding transformation of electromagnetic fields) can 
be mimicked by introducing electromagnetic materials with specific electromagnetic 
properties into the domain where the coordinates have been transformed [TJ [2j E]- As 
an example of such "coordinate-transforming" device an "invisibility cloak" has been 
suggested [HE]. In addition to perfect invisibility devices, similar approaches were 
applied to perfect lenses, description of the Aharonov-Bohm effect and artificial black 
holes. It has been known for a long time that the material relations of an isotropic 
magnetodielectric transform into certain bi-anisotropic relations if the medium is moving 
with a constant velocity (e.g., [5]). It was suggested that also this effect of time 
transformation can possibly be mimicked by a metamaterial [2], [6j [T] . 

We have recently introduced an alternative paradigm of transformation optics, 
where the required material properties are defined directly from the desired 
transformation of electromagnetic fields. This concept of field-transforming 
metamaterials was proposed in [5], where it was shown what metamaterial properties are 
required in order to perform the field transformation of the form E — > E = F(r, u;)Eq, 
Ho — > H = G(r, w)Ho, where F and G are arbitrary differentiable functions of position 
and frequency. In this concept we start directly from a desired transformation of 
electromagnetic fields and do not involve any space nor time transforms. A special case 
of similar transformations was proposed earlier as a numerical technique for termination 
of computational domain by modulating fields with a function decaying to zero at the 
termination [9]. The physical interpretation of this numerical technique in terms of a 
slab of a material with the moving- medium material relations was presented in [6] , and 
the pulse propagation in a slab of this medium backed by a boundary was studied in 

In this paper we generalize the concept of field-transforming metamaterials 
introducing methods to perform general bi-anisotropic field transformations, where the 
desired fields depend, in the general linear fashion, on both electric and magnetic fields 
of the original field distribution. Due to generality of the approach, the concept can be 
applied to a great variety of field transformations, not limited to cloaking devices. 
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2. Field transformations with metamaterials 



Let us assume that in a certain volume V of free space there exist electromagnetic fields 
E (r), H (r) created by sources located outside volume V (we work in the frequency 
domain, and these vectors are complex amplitudes of the fields) . The main idea is to fill 
volume V with a material in such a way that after the volume is filled, the original fields 
E and H would be transformed to other fields, according to the design goals. Here we 
consider metamaterials performing the general linear field transformation defined as 

E(r) = F(r, w )E„(r) + Ma(t, u)H (r) (1) 

V £ o 

H(r) = G(r, cu)H (r) + ,/^C(r, w)E (r). (2) 

V 

Here, scalar functions F(r,u), G(r,u) A(r,u) and C(r, oj) are arbitrary different iable 
functions. In this paper we consider only the case of scalar coefficients in the above 
relations, although the method can be extended to the most general linear relations 
between original and transformed fields by replacing scalar coefficients by arbitrary 
dyadics. 



2. 1 . Required constitutive parameters 

Substituting (JTJ) , (J2J) into the Maxwell equations and demanding that the original fields 
Eo and Ho satisfy the free-space Maxwell equations, one finds that the transformed 
fields E and H satisfy Maxwell equations V x E = — j'cjB, V x H = jujY) in a medium 
with the following material relations: 
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Relations (jSJ) and (jlj) describe a bi-anisotropic medium, whose constitutive relations can 
be conveniently written as [7] 

p • H + V e oMf + 3^) T ■ E, 
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Comparing with ([3]) and (jlj) we can identify the required material parameters of the 
field-transforming medium: 
The permittivity dyadic 

the permeability dyadic 

* = ^Jahc ( p2 + A2 ) 1+ ^Jg L ac^ a - AVF > x 1 < 8 > 

the nonreciprocity dyadic \ 

^=-FG^AC^ AG + CF ^ + 

+ ^IfO^AC l< AVC " FVG ' " < GVF " CV - 4 >' * 7 ' < 9 > 

and the reciprocal magneto-electric coupling dyadic (also called chirality dyadic) k 

1 1 1 



[-(AVC - FVG) - (GVF - CVA)] x 7. (10) 



^/e^2co FG-AC 
In the above equations J is the unit dyadic. 

3. Classification of generalized field-transforming metamaterials 

In this Section we classify the various cases of different materials that may need to be 
obtained depending on the required transformation (fTJ) , (j2]) . 

3.1. Reciprocity and nonreciprocity 

Considering the required permittivity and permeability we see that in the general case 
they both contain symmetric and anti-symmetric parts, so the material is nonreciprocal. 
The nonreciprocal parts of permittivity and permeability vanish, if coefficients A and C 
in (TjQ) and (j2J) equal zero, meaning that the desired field transformation does not involve 
magnetoelectric coupling, as in [8]. Note that in this case the medium can be still 
nonreciprocal due to its magnetoelectric properties, see below. The other important case 
when the permittivity and permeability are symmetric (actually scalar in the present 
case) is when the coefficients F and G vanish. These transformations define new electric 
field as a function of the original magnetic field and vice versa. Furthermore, we notice 
that if A = C = and, in addition, G = F, the wave impedance of the medium does 
not change. This is expected, because this is the case when the electric and magnetic 
fields are transformed in the same way. 

Next, let us consider what types of magnetoelectric coupling in field-transforming 
metamaterials are required for various transformations of fields. Here we will use the 
general classification of bi-anisotropic media according to [3 [TT] . This classification is 
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based on splitting the coupling dyadics x an d k> into three components, e.g. for k we 
write 

^= k7 + W + 1, (11) 

where A" and J denote the symmetric and antisymmetric parts, respectively. The 
scalar coefficient k equals to the trace of k, and it is the chirality parameter. In our 
case it is identically zero, meaning that field-transforming metamaterials (with scalar 
transformation coefficients) are nonchiral (their microstructure is mirror-symmetric). 
Also, the symmetric part A" is identically zero, and only the antisymmetric part 

K = FG^AC [ - (AVC ~ FVG) " (GVF " CV - 4)1 X 7 (12) 

remains. 



Table 1. Classification of reciprocal bi-anisotropic media. 



Coupling parameters 


Class 


K ^0, N = 0, 7=0 


Isotropic chiral medium 


K ^0, N^O, 7=0 


Anisotropic chiral medium 


K = 0, N^O, 7=0 


Pseudo chiral medium 


K = 0, N = 0, 7^0 


Omega medium 


K ^o, A 7 = 0, 7^0 


Chiral omega medium 


K = 0, N^0, 7^0 


Pseudo chiral omega medium 


K ^0, N^0, 7^0 


General reciprocal bi-anisotropic medium 



Classification of reciprocal bi-anisotropic media is given by Table 1 [7J. Field- 
transforming metamaterials possess reciprocal magneto-electric coupling as in omega 
media. The required effects can be realized by introducing inclusions in form of pairs of 
the letter Q, arranging them in pairs forming "hats" [TTI [T2] . 

Similarly, we split the nonreciprocity dyadic x as 

I = X 7 + ^ + £, (13) 

where Q and S denote the symmetric and antisymmetric parts, respectively, and x is 
the trace of dyadic x- From we can identify the parts of x as 

* = -re^ic ( - 4G + CF) ' (14) 

5=0, (15) 

5 = ^hTchic ^ AVC - FVG) - < GVF - GV4) i x 1 (16 > 

In contrast to the reciprocal magneto-electric coupling, in this case the trace of the 
dyadic x is generally non-zero. Looking at Table 2 [7J, we see that the transformation 
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Table 2. Classification of nonrcciprocal bi- anisotropic media. 



Coupling parameters 


Class 


X^O, Q = 0, S = 


Tellegen medium 


X^O, Q ^ 0, S = 


Anisotropic Tellegen medium 


x = 0,Q^0,S = 


pseudo Tellegen medium 


x = 0,Q = 0,S^0 


Moving medium 


x ^0,Q = 0,S^0 


Moving Tellegen medium 


x = 0,Q^0,S^0 


Moving pseudo Tellegen medium 




Nonreciprocal (nonchiral) medium 



medium can be either a Tellegen medium, a moving medium, or a moving Tellegen 
medium. 

If all the transformation coefficients do not depend on the spatial coordinates, 
the field-transforming metamaterial is the Tellegen medium [TT] , except the special 
case when AG = —CF. Tellegen media can be synthesized as composites containing 
magnetized ferrite inclusions coupled with small metal strips or wires p3] or as 
mechanically bound particles having permanent electric and magnetic moments [HJ[T5]. 



3.2. Losses and gain 

A bi-anisotropic medium is lossless if the permittivity and permeability dyadics are 
hermittian: 

f = I , V' — W (17) 
In addition, the magneto-electric dyadics for lossless media satisfy 

i+j*={i-m* (is) 

where * denotes complex conjugate [?]. 

Analysing formulas fl71)- f|T0l) for the required material parameters we see that if the 
transformation coefficients are complex numbers, in general the required material can 
be a lossy or gain-medium. For real-valued coefficients F, G, A, and C the permittivity, 
permeability, and reciprocal magnetoelectric coupling coefficient correspond to lossless 
media. However, the antisymmetric part of the nonreciprocity dyadic (describing effects 
of "moving" media) may correspond to lossy or gain media, depending on how the 
transformation coefficients depend on the position. 



3. 3. Presence or absence of sources 

In the derivation of the material relations in field-transforming metamaterials (jlj) 
we have assumed that in the medium there are no source currents, that is, the field 
equations are V x E = —ju)~B, V x H = joj~D. In addition, conditions 

V-D = 0, V-B = (19) 
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should be also satisfied in source-free media. For the general case we have no proof 
that conditions (fl9l) always hold. However, we have checked that they are satisfied 
in many important special cases, which include, for example, position-independent 
transformation coefficients and exponential dependence of the coefficients on the position 
vector. 



4. Particular cases 



Let us next consider some particular cases of field-transforming metamaterials to reveal 
the physical meaning of such transformations of fields. 

4-1. Backward-wave medium 

Consider the case of the field transformation E n — > H n — > ,/^-E, which 

corresponds to F = G = and A = C = 1 in ([IT),,(l2"]). Transforming electric field 
into magnetic field and vice versa, we reverse the propagation direction of plane waves 
in the transformation volume, so we expect that this would correspond to a backward- 
wave material filling the volume. And indeed we see that the material relations ([3]) and 
(BJ reduce to 

B = -fioK, (20) 
D = -e E, (21) 
which are the material relations of the Veselago medium [T6] . 



4-2. Field rotation 

Consider the case of the field transformation E n — > H n — > — which 

corresponds to F = G = and A — —1, C = 1 in ([I]),©. This kind of transformation 
relates to changing the polarization of the field: fields of a plane wave are rotated by 
—90 degrees around the propagation axis. The material relations and reduce to 

B = 0oH, (22) 

D = e E, (23) 

which are the material relations of free space. The reason why we do not see the change 
of polarization here is simply due to the fact that we wanted to consider F, G, A and 
C as scalars. The media (original or transformed) do not "feel" the polarization of the 
fields, so the rotation does not have any influence on the material relations. For instance, 
unpolarized light is obviously invariant under this transformation. We want to remind 
the reader that by considering the transformation coefficients F, G, A, C as dyadics, one 
can obtain more general transformations, which tailor also the field polarization in the 
general way. 
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4-3. Isotropic Tellegen medium 

Consider next a more general case when all the transformation coefficients do not depend 
on the position vector inside the transformation domain. In this case, gradients of these 
functions vanish, and we see that the material relations are that of an isotropic Tellegen 
material: 

D -»re^( G2 + C2 ) E -^§^H. (25) 

The medium is reciprocal only if x — ~ fg^ac = ^' ^ e can no * e ^ 
the transformation coefficients vary within the transformation volume, this generally 
requires simulation of a moving medium with the transformation metamaterial. 

4-4- Tellegen nihility 

Let us consider again the case where all the coefficients are position- independent, that 
is, 

VF = VG = VA = VC = (26) 

and also demand that the permittivity and permeability of the transforming medium 
become zero. This takes place when either 

A = jF, C = jG (27) 

or 

A = -jF, C = -jG (28) 
The material relations of the medium become 
B = -j^/e^E 

for the case (!27|) or 

D = Jv^oH 
for the case (128]) . 

These constitutive relations look similar to that of chiral nihility [17], but in this 
case the only non-zero material parameter is the Tellegen parameter. The Maxwell 
equations corresponding to relations (!30j) take the form 

V x H = -fc H, V x E = fc E, (31) 



Generalized field-transforming metamaterials 



9 



where ko = uj-Je^Iq is the free-space wavenumber. Thus, solutions of Eq. (13TI) are 
eigenfunctions of operator rot=V x /. If d/dx = d/dy = then for the waves 
propagating in ^-direction 

E x = -e exp (-jk z) H x = e^n' 1 exp {-jk z) 

E y = je exp (-jk z) H y = je ^ _1 exp (-jk z) (32) 
E z = H x = 

where 77 is the wave impedance. Substituting expressions (I27j) or (12"51) for A and C into 
the original definitions (pQ),(j2]), we obtain the following formula for the wave impedance: 
F 

V = ^Vo (33) 

where 770 is the wave impedance of vacuum. 

As in the case of chiral nihility [T7] and the Beltrami fields in chiral media [HJ HE] , 
the fields f)32p correspond to circularly-polarized waves. However, in contrast to those 
two cases, where the phase shift between electric and magnetic fields equals exactly to 
7r/2, in the considered case it depends on relation between F and G according to (J3"3"j) . 
If F = G, the electric and magnetic fields oscillate in phase. 

4-5. Position-independent transformation coefficients 

An interesting general conclusion can be drawn with regard to field transformations 
with all the coefficients being position-independent. In this case the wavenumber in the 
field-transforming metamaterial is the same as in the original medium (free space in our 
case). This can be shown by using the equation for the wavenumber A; in a Tellegen 
medium fTj 



k = k sje r /j, r - x 2 , (34) 

where e r and fi r are the relative permittivity and permeability. Since 

2 _ (G 2 + C 2 )(F 2 + A 2 ) (AG + CF) 2 _ 
£r/ir X (FG - AC) 2 (FG-AC) 2 ' 1 j 

( 1341 simplifies to 

k = ko. (36) 

This means that for transformations of this class only the wave impedance of the 
transforming medium needs to be changed, but not the wavenumber. One can notice 
that this property indeed holds in all of the above special cases where the coefficients 
are position-independent. 

4-6. Moving omega medium 

Next we consider another class of field-transforming media assuming A(r,u) = 
C(r,u) = 0. Transformations of this type correspond to moving omega materials with 
the constitutive relations of the form: 

D = £o fE-^VGxH, B = /i gH + ^VFxE. (37) 
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Both the nonreciprocity dyadic x an d the reciprocal magnetoelectric coupling dyadic k 
are antisymmetric: 



j [VF VG\ = _ 1 /VF VG\ = 

x = ~k VT + -g) x j ' k = ~2k [-T ~ -g) x 7 ' (38) 

where the first quantity describes moving-media effects and the second one corresponds 
to omega-medium properties. If the two transformation coefficients are equal (F = G), 
the omega coupling vanishes. 

The Maxwell equation for the transformed fields read: 



V x E(r) = -j^ ggH(r) + ^VF(r) x E(r) 

V x H(r) = j^oggE(r) + ^VG(r) x H(r). 



(39) 



4-7. Artificial moving medium 



Let us consider a more special case when the transformation preserves the wave 
impedance of the medium [F(r) = G(r)] [8]. Material relations (137|) take the form 
of material relations of slowly moving media (velocity « < c, where c is the speed of 
light in vacuum), see, e.g., [19J: 

D = -j^— (e E+ x H) « e E+ x H 



■sjl-v 2 /c 2 

B = (moH -1vxE)^„H-1vxE 



(40) 



In the last relation we have neglected the second-order terms v 2 /c 2 . The "velocity" of 
the transforming medium can be identified by comparing with (|37|) as 

V = "J r- 41 

oj F 

To get more insight, let us analyze the one-dimensional case considering 
transformation functions depending only on one coordinate z and assuming that 
F{z) = G(z) is an exponential function of z. Two cases should be distinguished: 

• VF/F is imaginary, e.g., F(z) = e~i az with a real parameter a. In this case the 
effective "velocity" is real. This metamaterial simulates moving media. Note that 
the required material relations have the form of the relations for slowly moving 
media even if the equivalent velocity is not small. 

• F(z) is a purely real function, e.g., F(z) = e~ az , where a is real. In this case the 
effective "velocity" of the medium is imaginary. Note that the material relations of 
a medium traveling faster than light formally have an imaginary vector coefficient in 
the second term of fj4*0l . However, one cannot say that this metamaterial simulates 
media moving faster than light, because in that case also the effective permittivity 
and permeability would be imaginary quantities. 

Let us examine the source-free conditions V • D = and V • B = for this special 
case of F = G = e~ az . Actually, = = ~ a z °- Then [t follows from ^> 

V • D = V • [F(r)e E (r)] + • [az x F(r)H ] (42) 

u 
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Simple calculation using the free-space Maxwell equations relating Eo and Ho shows 
that V ■ D = 0. Condition V ■ B = can be checked similarly. 

Let us check the passivity condition for this medium. The period-average of the 
time derivative of the field energy is expressed as [7] : 

ldW\ joo i 



dt L 4 



■f* ■ (M - M f ) • f (43) 



where 



i=( E M = ( 601 _ v 7 ^^ 

y^ojM)X Vol 




(44) 



-e~ 2az > (45) 
V 



For real transformation coefficients we take as an example F(z) = e az , where a > 
For a plane wave propagating along z, calculation leads to 

ldW\ 

\~dT/ t 

This implies that the medium possesses losses. However, for a plane wave traveling in the 
opposite direction (along —z), the sign in the above expression is reversed, corresponding 
to a gain medium. This is expected because in the last case the field amplitude grows 
along the propagation direction. In the case of imaginary a this value is equal to zero 
and the medium is lossless. 



5. Transmission through a slab of a field-transforming medium 

Let us consider an infinite slab of a field-transforming metamaterial in vacuum, which 
is excited by a normally incident plane wave. Here we continue studying the case of 
moving-medium transformations, that is, F = G and A = C = 0. Let the thickness 
of the slab be d, and the field-transforming properties be defined by a function of one 
variable z — the coordinate in the direction normal to the interfaces: F(z) = G(z). 
The general solution for the eigenwaves in the slab is, obviously, 

E(z) = F(z) (ae- jkoZ + be jh ° z ) (46) 

where E is the electric field amplitude, and a and b are constants, because this is how 
the free-space solution is transformed in this metamaterial slab. If F = G, the wave 
impedance in the transforming medium does not change. 

Writing the boundary conditions on two interfaces, reflection and transmission 
coefficients can be found in the usual way. As a result, we have found that the 
reflection coefficient does not depend on function F(z) at all, and it is given by the 
standard formula for a reflection coefficient from an isotropic slab at normal incidence. 
In particular, if the permittivity and permeability of the slab equal to that of free space, 
the reflection coefficient R = 0. The transmission coefficient in this matched case reads 
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where F(d) and F(0) are the values of the transformation function at the two sides of 
the slab. 

The non-reciprocal nature of field-transforming metamaterials is clearly visible here. 
If function F(z) is real and it decays with z, then fields of waves traveling in the positive 
direction of z decay, while fields of waves traveling in the opposite directions grow. 
Obviously, such performance cannot be realized by passive media (see the calculation of 
the absorbed power in Section H~T1) . If F(d) = F(0) = 1, the slab is "invisible", although 
inside the slab the field distribution can be pretty arbitrary, as dictated by the given 
function F(z). 

6. Conclusions 

Two different approaches to design of metamaterials which transform fields in a desired 
way are known. The first one has been proposed and developed in [fl [2], and it is 
based on coordinate transformations. The second one has been introduced in [8] and 
developed here, and it transforms the fields directly, according to a desired prescription. 
This technique can be potentially applied to a very wide variety of engineering problems, 
where a certain field transformation is required. One of the most important property of 
the field transformation approach is that a medium with prescribed parameters turns 
out to be nonreciprocal with except of some special cases. Most often, a metamaterial 
simulating moving media is required. Potential physical realizations of metamaterials 
with the properties required for general linear field transformations (such as artificial 
moving media) are discussed in [HI [3, [EU EE ED] . 
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